arXiv:1504.02820vl [astro-ph.GA] 11 Apr 2015 


Mon. Not. R. Astron. Soc. 000, [T| fl4| (2007) Printed 14 April 2015 (MN IATeX style file v2.2) 


The soft and hard X-rays thermal emission from star 
cluster winds with a supernova explosion 


A. Castellanos-Ramirez, 1 * A. Rodriguez-Gonzalez , 1 A. Esquivel , 1 
J.C. Toledo-Roy 1 , J. Olivares , 2 & P.F. Velazquez 1 * * 

1 Instituto de Giencias Nucleares, Universidad Nacional Autonoma de Mexico, Apdo. Postal 7 0-543, 04510, Mexico, D.F., Mexico 
2 Univ. Grenoble Alpes, 1PAG, F-38000 Grenoble, France CNRS, IPAG, F-38000 Grenoble, France 


Draft Version, 14 April 2015 


ABSTRACT 

Massive young star clusters contain dozens or hundreds of massive stars that inject 
mechanical energy in the form of winds and supernova explosions, producing an out¬ 
flow which expands into their surrounding medium, shocking it and forming structures 
called superbubbles. The regions of shocked material can have temperatures in excess 
of 10 6 K, and emit mainly in thermal X-rays (soft and hard). This X-ray emission is 
strongly affected by the action of thermal conduction, as well as by the metallicity of 
the material injected by the massive stars. We present three-dimensional numerical 
simulations exploring these two effects, metallicity of the stellar winds and supernova 
explosions, as well as thermal conduction. 

Key words: Hydrodynamics - galaxies: star clusters: general - (Galaxy:) open clus¬ 
ters and associations: general - ISM: bubbles - shock waves - stars: massive - stars: 
winds, outflows - X-rays: ISM 


1 INTRODUCTION 


It is well known that masive O-B type stars inject a con¬ 
siderable amount of mechanical energy into the interstellar 
medium (ISM), in form of stellar winds or supernova (SN) 
explosions. The energy input by these events is sufficient to 
drive strong shocks that expand into the ISM generating a 
structure called bubble. 


The model proposed by Weaver et al. (19771 and later 


expanded by Chu & Mac Low (19901 and Chu et al. (19951, 


is considered the standard model of bubbles driven by stellar 
winds. It considers the injection of mechanical energy to the 
ISM from stellar winds that results in the formation of a 
bubble. This bubble is surrounded by a cool shell of ISM 
material that has been swept by the expanding shock front. 
The shocked (and thereby heated and compressed) material 
in the interior of the bubble emits considerably in X-rays, 
whereas the outer, cooler shell emits at optical wavelenghts. 


The original Weaver et al. (1977) model considers a sin¬ 


gle stellar wind source. Some time later, in order to explain 
what is now known as superbubbles, the model was extended 
to include multiple wind sources (see|Chu fe Mac Low|199~0) 


Chu et ah|l995 Canto et al.|2000 Silich et al.|2004l . 

The simplest model of superbubble formation is as fol¬ 
lows. Consider a cluster with N stars each having different 


* antonio.castellanos@nucleares.unam.mx 


mass-loss rate M Wi i and a wind velocity u w ,i. Since the stars 
inject mechanical energy in the form of stellar winds, the 
total mechanical luminosity is given by 


Iw 



i =1 


At first, the stellar winds collide with each other and 
with the environment inside the cluster radius. Thus the 
space between the stars is filled with hot shocked material 
from the winds. This happens until a stationary flow is es¬ 
tablished, giving rise to a common cluster wind that forms a 
super shell. As this supershcll expands through the surround¬ 
ing ISM it creates a superbubble structure with the follow¬ 
ing structure (Weaver et al. 1977 Rodrfguez-Gonzalez et al. 


2011 Velazquez et al.||2013|): 


(i) The innermost region located near the stars (where 

their winds collide) produces thermal hard X-ray emission 

(if the stellar winds have terminal velocities larger than 1000 
km s -1 ), and driving the expansion of the bubble through a 
pressure difference between the hot and dense interior and 
the colder and less dense environment. 

(ii) After the individual winds from the stars coalesce into 

a cluster wind, it expands freely from the cluster radius out¬ 
wards. In this zone X-ray emission is important only close 
to the cluster radius, and it consists mostly of soft X-rays. 

(iii) Behind the main shock pushing into the ISM a re- 
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2 Castellanos-Ramirez et al. 


verse shock is formed. The reverse shock encounters the 
freely expanding wind and compress and heats it to soft 
X-ray emitting temperatures. The region filled by shocked 
wind is quite extended and dominates the emission in X- 
rays, particularly in the soft energy bands. 

(iv) The outermost region of the supperbubble consists 
of a shell of shocked ISM that has been swept up by the 
main shock. Beyond this zone there is only unperturbed ISM 
material. 


The original wind blown bubble (WWB) model pro¬ 
posed by Weaver et al. (1977) overpredicts the X-ray lumi¬ 
nosity. One reason is that this model includes thermal con¬ 
duction and in consequence produces a denser interior that 
in the case without it, which in turn increases the X-ray 
luminosity. Furthermore, the wind blown bubble models do 
not take into account the radiative losses within the clus¬ 
ter radius, and this can have a significant impact on the 
luminosity (see Rodriguez Gonzalez et al. 2011). Recent ex¬ 
amples of this are the works of Dunne et al. (2003) and 
Reyes-Iturbide et al. (2009), which predict X-ray luminosi¬ 
ties that exceed that of the observations by about one order 
of magnitude. 

On the other hand, there are others models that predict 
an X-ray emission that underestimates the observed values 
(for instance, see the work of Harper-Clark & Murray 2009; 
Rogers & Pittard 2014, in which only the cluster wind region 
is considered), a problem for which different solutions have 
been explored. |Chu fe Mac Low (1990I proposed that in or¬ 
der to increase the luminosity of X-rays (so as to match the 
observations) one should consider shock waves produced by 
the explosion of supernovae inside the star cluster. |t evens | 


& Hartwell (2003) presented models where the luminosity of 


soft X-rays is obtained as a function of the mass loss rate, 
the cluster radius and the wind terminal velocity. They do 
not take into account mass loading, but they consider it 
can be relevant for the study of soft X-rays in this type 
of massive clusters. The work of Silich et al. (2001) deals 
with the effect on the X-ray emission of the high metal con¬ 
tent injected by the massive stellar winds and the SN rem¬ 
nants. Rodrfguez-Gonzalez et al. (2011) showed that super¬ 


novae occurring near the centre of the cluster are not capable 
of reproducing (completely) the luminosity observed in X- 
rays, and neither do they help explain the kinematics of the 
shell (without consider thermal conduction). They instead 
showed that off-centre SN explosions (for N70 and N185, see 
also Reyes-Iturbide et al. 20141 could help explain the two or 
three orders of magnitude difference between the luminosity 
observed and the standard model predictions. However, in 
these models the X-ray luminosity agrees with the observed 
value for only ~ 10 000 years, making the probability of ob¬ 
serving them in this regime rather low. On the other hand, 


Velazquez et al. (2013) presented models of the M 17 super¬ 


bubble where they considered the contribution of the gas of 
the parental cloud in the evolution. They showed that the 
mass loading from the parental cloud can help increase the 
luminosity of soft X-rays by up to an order of magnitude. 
In such work neither the metallicity nor thermal conduction 
were considered. 

The Large Magellanic Cloud (LMC) is filled with su¬ 
perbubbles with important soft X-rays emission. Some of 
these superbubbles (for instance, DEM L50 and DEM L152, 


Jaskot et al. 2011) show evidence for off-centre super¬ 


nova events which seem to interact with the external shell 
pushed by the stellar winds. In the observations made by 


Jaskot et al. (2011) the gas of the supernova remnant is still 


seen. This remnant is located close to the superbubble edge. 
These objects have luminosities up to an order of magnitude 
higher than those predicted by the model of |Weaver et al. | 
(19771. 

Moreover, the numerical models that appear in | Jaskot] 
et al. (2011) produce luminosities that are two orders of 


10 36 erg s 1,1 


magnitude lower than the observations (^ 

Therefore the authors explored the effect of the metallicity 
and of mass loading by clouds to bridge the luminosity deficit 
in soft X-rays. They calculated the mass of metal enriched 
material injected by the supernova explosions ( Maeder|1992 


Oey 1995 Silich et al. 2001 Anorve-Zeferino et al. |2009 


and found that metallicities from 3 to 10 times solar can 
be achieved, and using the equations of Silich et al. (2001), 
they concluded that the effect is not sufficient to account 
for the differences. They conclude that the main mechanism 
that can explain such an important enhancement of the total 
X-rays luminosity is mass loading. 


Recently, Rogers & Pittard (2014) presented a study 
of the soft X-ray emission during the various evolutionary 
stages of massive stars embedded in a dense giant molecular 
cloud (GMC), going through the red supergiant and Wolf- 
Rayet stages up to the supernova phase. They showed that 
the inclusion of the GMC results in a short lived attenuation 
of the X-ray emission of the cluster, during the time before 
an important fraction of the material is carried away from 
the wind interaction region. After this occurs, the luminosity 
remains practically constant. 

The X-ray emission of a star changes substantially as it 
goes through distinct evolutionary stages. For instance, the 
X-ray luminosity drops abruptly during the red giant phase 
and increases substantially once in the Wolf-Rayet phase. 


Rogers & Pittard (20141 show that, in spite of the differ¬ 


ences between their models and some observations, their re¬ 
sults agree reasonably with other observations, such as the 
case of M17 and the Rosette Nebula. They found that the 
emission produced by their model during the early wind- 
dominated phase is smaller compared to the prediction from 


the standard model (Weaver et al. 1977 Chu & Mac Low 


19901, but larger than the emission expected in models that 


only consider the emission at the interaction region of the 
winds of massive stars (the cluster wind). Finally, for stars 
in the main sequence, they found luminosities two or three 
orders of magnitude above those predicted by the standard 
model, lasting for more than 4.5 kyr. 


In this work, we present a series of numerical models 
exploring the effects of supernova explosions, metallicty and 
heat conduction in the thermal soft and hard X-rays lumi¬ 
nosity of a massive star cluster. The paper is organised as 
follows: in Section 2 we present the numerical setup of our 
models, describe the implementation of the thermal conduc¬ 
tion and the metallicity in the gas dynamics equations, and 
in Section 3 we show the resulting synthetic emission in the 
soft and hard X-ray bands as well as a brief discussion of 
our results . In Section 4 we have made some comparisons 
of our numerical models with four interesting observed bub¬ 
bles. Finally, a summary is given in Section 5. 
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X-ray thermal emission in cluster winds with a SN 3 


2 THE NUMERICAL SIMULATIONS 


With the purpose of exploring the effects of the interaction 
and influence of the SN explosions and metallicity, as well 
as thermal conduction in the cluster stellar winds, we per¬ 
formed a series of numerical simulations, and estimated the 
soft and hard X-ray emission that would be produced. 

We used the huacho code (see Esquivel et al.||2009 & 
Raga et al.|2009 1 to perform all numerical simulations. The 


code solves the hydrodynamic equations ( 1)3 1 on a three 
dimensional uniform Cartesian mesh, using a second order 
finite volume method with HLLC fluxes (Toro et al. 1994) 
and a piecewise linear reconstruction of the variables at the 
cell interfaces with a minmod slope limiter. The code also 
includes radiative losses and isotropic thermal conduction: 


g+V-(pu) = 0, 


+ V ■ (puu + IP) = 0, 


djpu) 
dt 

^ + V ■ [u (E + P)] = L rad (Z, T) + V ■ q, 


(1) 


( 2 ) 


(3) 


where p, u, T, P and E are the mass density, velocity, tem¬ 
perature, thermal pressure and energy density, respectively, 
I identity matrix, 7 is the heat capacity ratio, L rad (-Z, T) 
is the energy loss rate, and q is the heat flux due to elec¬ 
tron conduction (see sub-section 2.2). The system is closed 
with an ideal gas law given by E = p|u| 2 /2 + P /(7 — 1). To 
find the energy loss rate, we use a tabulated cooling func¬ 


tion from the freely available Chianti database (Dere et al. 
|1997| |Landi et al.||1996| ) . As we show in Figure 0 we have 
constructed a cooling function with a range of metallicities 
of 0.1-30 Zq. 

The computational domain is a cube of 140 pc on a side, 
discretised by 256 3 cells in a uniform grid, yielding a reso¬ 
lution of 0.5469 pcQ From the number of massive stars, one 
can estimate the mass of the star cluster (~ 3500 M@ using 
Starburst99, Leitherer & Heckman (1995). In our model we 
did not consider the total mass of the paternal cloud, but 
we selected the size of the simulation box so that it contains 
a typical superbubble of radius (~50 pc, i.e. DEM L 50 and 
DEM L 152). 


The simulations include 15 stellar wind sources placed 
randomly within the cluster radius (R c = 10 pc, and the 
distribution is the same for all the simulations). The stellar 
winds are imposed in spherical regions of radius P w = 6.03 x 
10 18 cm (1.95 pc), corresponding to 5 pixels of the grid, and 
have a temperature T w = 10° K. All the stars have the same 
mass-loss rate, M w = 10 ~ 6 M@ yr -1 , and wind velocity of 
u w = 1500 km s _1 . We turned on the stellar winds at the 
beginning of the simulation. The rest of the computational 
domain was initially filled by a homogeneous enviromnent 
with temperature To = 10 4 K and density no = 2 cm -3 . 

We impose a SN inside the bubble formed by the winds 
at four different times: 2, 3, 5 and 7.5 x 10 5 yr (each cor¬ 
responding to a different model). These times were chosen 


1 We have tested that the resolution is sufficient in terms of the 
numerical convergence of the overall features of the simulations. 
However some details of the flow and the exact values of the 
luminosities still depend on the resolution; see the discussion at 
the end of Section 3. 



T [I ] 


Figure 1 . Cooling function for a range of metallicities between 
0.1 to 30 Zg. 


to control the distance from the site of the supernova to the 
supperbuble shell. Since we do not follow the evolution of 


the star cluster that produces the shell (as Rogers & Pittard 


20141, these times are related to the superbubble dynamical 


age that would be observed and not with the star cluster 
age. The superbubble dynamical age is smaller than that 
of the stars because it does not include the time needed to 
clear up the material between the stars and form a com¬ 
mon bubble. While in most of the models the supernova is 
placed at the centre of the star cluster, we have included two 
off-centre models: one in which the SN is placed 5 pc from 
the centre, and one where it is near the edge of the bubble 
(10 pc from the centre). The supernova explosion was im¬ 
posed by injecting a total energy of lxlO 51 erg and 2 Mg 
of mass in a region with a radius of 2 pc. Half of this energy 
was injected as kinetic energy (with velocity following an 
increasing linear profile with radius, and constant density 
and temperature inside the imposition region), and the rest 
is thermal energy ( Toledo-Roy et al.|2014 l. 

To explore the effect of metallicity we have performed 
some runs with a homogeneous metallicity for all the three 
components (the ISM, the stellar winds and the SN) and 
some models with a different metallicity for each of these 
components. In the homogeneous metallicity models we have 
used a metallicity of 0.3 Zq. For the variable metallicity 
models, following |Silich et al.pOOT) we use Zism = 0.3 Zq 
for the ISM, Z„i n d = 3.0 Zq for the mass injected in the 
form of winds, and Zsn = 10 Zq for the SN ejecta. We have 
also included thermal conduction in two of the models. 

The parameters of the simulations are listed in Table [Tj 
As can be seen from the table, we named the models to re¬ 
flect the parameters used: the number after ‘sr’ corresponds 
to the location of the SN in pc; it is followed by ‘tsn’ and 
another number to indicate the time of the SN detonation 
since the winds sources were turned on (in units of 10 s yr); 
next there is the letter ‘z’ followed by the number 0.3 for the 
uniform metallicity models or the letter ‘v’ for the variable 
metallicity runs; for the models with thermal conduction a 
letter ‘C’ is appended at the end of the model name. 
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Table 1 . Parameters of the simulations. 


Model 

SN locus 

[pc] 

SN detonation Metallicity 
time 

[10 s yr] Zq 

Thermal 

Conduction 

sr0tsn2z0.3 

0 

2 

0.3 

no 

sr0tsn3z0.3 

0 

3 

0.3 

no 

sr0tsn5z0.3 

0 

5 

0.3 

no 

sr0tsn2zv 

0 

2 

variable 0. 

no 

sr0tsn3zv 

0 

3 

variable 0, 

no 

sr0tsn5zv 

0 

5 

variable 0, 

no 

sr0tsn7.5zv 

0 

7.5 

variable 0, 

no 

sr0tsn2zvC 

0 

2 

variable 0, 

yes 

sr0tsn5zvC 

0 

5 

variable 0, 

yes 

sr5tsn5zv 

5 

5 

variable 0, 

no 

srl0tsn5zv 

10 

5 

variable 0, 

no 


a For these models the metallicity of the ISM is 0.3 ZQ, 3 Zq 
for the stellar winds, and 10 Zq for the SN ejecta. 


2.1 Adding the effect of metallicity to the cooling 


The cooling in the code is added as a source term after 
updating the hydrodynamic variables. At the end of each 
timestep, we estimate the cooling by interpolating a tabu¬ 
lated cooling curve which is, for a given metallicity, a func¬ 
tion of the temperature. The energy loss is then subtracted 
to the internal energy of each cell at every timestep. 

For the runs with a uniform metallicity this is a simple 
linear interpolation (in temperature) of a single table that is 
generated by the CHIANTI database. For the runs with vary¬ 
ing metallicity we created a series of tables for metallicities 
in the range of 0.1-30 Zq; these are plotted in Figure [l] 
Along with the gas-dynamic equations (eqs. |TJ|3j) we con¬ 
sider the metallicity Z as a passive scalar by including an 
extra equation of the form 


dZp 

dt 


+ V ■ (Zpu) = 0. 


(4) 


Using the metallicity value at each cell we do a bi-linear 
interpolation (with metallicity spaced linearly and temper¬ 
ature logarithmically) to estimate the cooling to be ap¬ 
plied there. For the ISM the metallicity is set to 0.3 Zq 
at the start of the simulation. For the winds and super¬ 
novae the gas is injected into the simulation either with 
•Zwind = Zsn = 0.3 Zq (for the homogeneous models) or 
Z w ind = 3 Zq and Zsn = 10 Zq for the inhomogeneous 


models. 

The average metallicity in each region can be calculated 
(using Silich et al.|2001) as: 


Z = 


Af z ,ej + AT,, ism 

M e j + M; sm 


(5) 


where, M Zje j and A/ Z ,i sm are the masses of the metallic ejecta 
(by winds and/or SN) and swept up interstellar gas, respec¬ 
tively, while M e j and Mi am are the total masses of the ejected 
and swept up interstellar gas. 

In general, the most metal enriched regions are found 
behind the contact discontinuity that separates the main 
and reverse shocks. Even though some mixing occurs at the 
interface (mainly due to hydrodynamical instabilities and/or 
turbulence), since the swept up ISM mass is larger than the 


ejected mass the metallicity of the shell remains close to that 
of the ISM. 


2.2 Thermal conduction 

In order to include the effect of thermal conduction by free 
electrons in our numerical simulation, we add a heat flux 
term (V • q) in the right hand side of the energy equation 

<§. 

The heat conduction due to collisions with free electrons 
in a plasma is given by the classical |Spitzer : (| 1962[ | law: 

q = - fcVT, (6) 

where k is the thermal conductivity given by 

k = /3 T s/2 . (7) 


where, for a fully ionized hydrogen plasma (3 « 6 x 


10 


7 erg s 1 K cm 1 (see Spitzer 1962). The result relies on 


the assumption that the mean free path is small compared 
to the scale-length of temperature variations (A -C T /j VTj). 

When the mean free path of the electrons is compa¬ 
rable or larger than temperature scale-length the heat flux 
saturates. In this regime the heat flux can be estimated by 
the local sound speed (c s ) and pressure ( P ), as described by 
Cowie & McKee ( 1997| ): 


Qsat — 5<f> sCs P, (8) 

where (j> s is a factor of order unity (we have used (p s = 1.1). 

At every timestep we compute the heat fluxes in the 
classical and the saturated regimes, keep the smaller one 
and introduce its divergence as a source term to the energy 
equation. We have to mention that the thermal conduction 
timescale is smaller than the hydrodynamic one determined 
by the standard CFL condition. For this reason we apply a 
sub-stepping method to include the source term (we take on 
the order of 100 sub-steps to integrate the source term for 
each hydrodynamical step). 
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Figure 2. Thermal soft (0.1—2 keV, in red lines) and hard (2- 
10 keV, in blue lines) X-ray emission coefficients for a range of 
metallicities between 0.1 to 30 Zq. 


2.3 X-ray emission coefficients 


We take the output from the hydrodynamical simulations 
to estimate the X-ray luminosity in all the models. We con¬ 
sider that the emission coefficient in the low-density regime 
is j v (n, Z,T) = rieX(Z,T), where n e is the electron density, 
and X {Z, T) is a function of the the temperature (T) and the 
metallicity (Z). For a given metallicity, the function x can be 
computed and integrated over an energy band using the chi- 
anti atomic database and its associated IDL software ( |Dere| 
et al.|1997 Landi et al.|1996 |. We have computed x(Z, T) for 


various metallicities (Z = 0.1, 0.3, 1, 3, and 10 Z@), using 
the ionisation equilibrium model by |Mazzota et ah] (1998|), 
over a range of temperatures from 10 4 to 10 y K. The emis¬ 
sion coefficients were integrated over two energy bands: soft 
X-rays (0.1-2 keV), and hard X-rays (2-10 keV). The result 
is a two dimensional table of coefficients that is function of 
temperature and metallicity. Figure [2]show the thermal soft 
(red lines) and hard (blue lines) X-ray emission coefficients, 
respectively, as function of temperature for several metallic¬ 
ities. 


From the results of the simulations we obtain the den¬ 
sity, temperature and metallicity in every computational cell 
and perform a bilinear interpolation to get y, and then use 
it to determine the emissivity for the two energy bands. The 
contribution of all cells are then added to compute the total 
X-ray luminosity both for soft and hard X-rays. 


3 RESULTS 

The colour maps of Figures [3] and [4] show the density and 
temperature at three different evolutionary times. These 
were chosen to show the effect of the SN explosion in the 
X-ray emission (see the next section). We present a time 
slightly before the SN explosion (top row), at the peak of 
luminosity after the explosion (middle row) and once the to¬ 
tal luminosity has diminished back to a value near its pre-SN 
level (bottom row). The columns correspond to three differ¬ 
ent models: in the left and central columns the SN ocurrs at 
the centre of the cluster, without thermal conduction, and 
with thermal conduction respectively, and in the rightmost 


column the SN is 10 pc off-centre (the position of the SN 
is indicated by a star in the top row). Following the time 
sequence in the columns of this figure, it can be seen that 
the SN ejecta reach the edge of the wind bubble and push 
it further into the ambient medium. Due to the particular 
position of the stars in these models, the gas distribution 
inside the wind bubble favors the expansion of the SN eject 
towards the upper right corner of the simulation box, and 
thus the blowout is more pronounced in this direction, the 
effect being larger if the SN explodes off-centre (at the edge 
of the star cluster; see the rightmost panels). 


3.1 Soft X-ray emission 


Following the pressure driven model discussed by |Chu &| 
[M ac Low]( 1990|, the soft X-ray luminosity can be estimated 
from 


L x = 3.29/(t)£Z/37 
where 

U \ 125 = 1/2 

^ T = ~33 5r ^ 

with 


33/35^17/35^19/35 


n. 


0 


[erg s" 


, 5 3 

H- T - T 

3 11 


11/3 


r = 0.16 L 3 - 8 / 35 n 0 - 2 / 35 t® /35 , 


(9) 


( 10 ) 


( 11 ) 


£ is the gas metallicity, L 37 = L w /10 3 ' where L w is the 
mechanical luminosity of the cluster (in erg s _ ), no the 
interstellar medium density and te is the cluster lifetime 
in Myr. In all the models presented here, the mechanical 
energy injected by the winds is 1.1 x 10 3 ' erg s -1 . With this 
mechanical energy the total X-ray luminosity for the s tellar 
10 33 erg s _1 (see eq.pll and also Chu 


wind contribution is ■ 


& Mac Low 1990), for an interstellar medium density of 2 
cm ' 3 and metallicity of 0.3 Zq after an evolution time of 
2 x 10 s yrs. 

We have computed the soft X-ray emission for all the 
models at 10 4 yr intervals. Figure |5| shows the evolution 
of the total soft X-ray luminosities for all models without 
thermal conduction and with the SN placed at the centre of 
the star cluster. The red lines show the models with uniform 
metallicity, and the blue ones with variable metallicity. 

The X-ray luminosity before the supernova event is in 
agreement with the value predicted by |Chu fe Mac Low| 
(1990). Shortly after the supernova explosion the luminos¬ 


ity increases dramatically. We calculated the time interval 
in which the soft X-ray luminosity remains above 10 34 , 10 35 
and 10 36 erg s _1 (Af s , 34 , At s ,35 and Af s , 36 , respectively). 
The maximum luminosity achieved and these time intervals 
are shown in Table [2] 

From Figure [5] we can see that the general shape of the 
luminosity curve after the SN explosion is quite similar in 
all the models. All models reached the same maximum lu- 


3 x 10 35 erg" 

a4 


and have At s ,34 ~ 6 x 10 4 yr 


minosity of - 

and At Sj 35 ~3xl0 4 yr. In these models, in which the SN 
explosion occurs at the centre of the stellar distribution, lu¬ 
minosities above 10 36 erg s _1 are never reached. Our results 
show only small differences in the soft X-ray emission be¬ 
tween models with uniform metallicity and those with vari¬ 
able metallicity. The similarity of the emission across the 
models indicates that the emission is dominated by swept 
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Figure 3. Density maps at different stages of the evolution and for different models. The columns correspond to three distinct models: 
in the panels of the left column, the SN explosion occurs at the centre of the cluster (model sr0tsn5zv; panels a, d and g); those of the 
central column show the same model but including thermal conduction (sr0tsn5zvC; panels b, e and h); and in those of the right column 
the SN explosion occurs 10 pc off-centre (srl0tsn5zv; panels c, f and i). For all models shown here, the SN explosion occurs at t = 5 X 10 5 
yr and the metallicity of the gas varies across components, as discussed in the text. The rows correspond to three relevant evolutionary 
stages: just before the SN explosion (top row), at the peak of X-ray luminosity (middle row), and once the luminosity has approximately 
returned to its original level (bottom row). The position of the SN is marked with a star in the panels of the top row. The spatial scale 
is the same across all panels and is shown in panel g. 


up ISM material. This is partly because the high metallic¬ 
ity gas (winds and SN) is kept at a temperature too high 
for thermal soft X-rays to be important. Thermal conduc¬ 
tion allows energy transport with reduced bulk motions, re¬ 
sulting in a denser inner region. Weaver et al. (1977) esti¬ 
mated a significant increase in the soft X-ray luminosity (up 
to 2 orders of magnitude) with respect to models without 
thermal conduction. Figure [6] displays the evolution of the 


soft X-rays luminosities for the 2 models with thermal con¬ 
duction (sr0t2e5zvC and sr0t5e5zvC, magenta lines), and 
their counterpart without thermal conduction (sr0t2e5zv 
and sr0t5e5zv, blue lines). As can be seen, thermal conduc¬ 
tion does increase the maximum luminosity of the models, 
but only by a factor of ~ 1.25, and the luminosity returns to 
a value a factor of 2 larger after the SN explosion. The time 
of emission above 10‘ u erg s -1 also increases by a similar 
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Figure 4. Temperature maps, corresponding to the same panels as in Figure [3] 


factors of 1.2 and 1.6 for the SN imposed after 2 x 10 5 yr 
and 5 x 10 5 yr, respectively. Approximately the same time- 
span increase is found for emission above 10 3 ° erg s^ 1 ; see 
Table[2] These increments in the time interval with emission 
above 10 34 and/or 10 35 erg s _1 will enhance the chance of 
such luminosities being observed. 

We can see that the inclusion of different metallicities 
and/or thermal conduction induces only small discrepancies 
in the soft X-ray emission. 

From this models it is clear that the supernova explo¬ 
sions are a crucial ingredient for the thermal X-ray emission. 
The presence of a SN can explain the extra X-ray luminosity 
observed in several superbubbles. However, when the super¬ 
nova event occurs in the centre of the cluster, the soft X-ray 
luminosity only reaches a few times 10 35 erg s -1 , still falling 


short of some of the observed values (e.g., those of Jaskot 
et al.| 2011). We find that models with a centred explosion 


seem to still be underluminous. 

A possibility that results in luminosities above 
10 36 erg s _1 is to place the SN at a distance from the centre 
of the star cluster. For these reason we have included models 
sr5tsn5zv and srl0tsn5zv where the SN occurs at 7 ?sn = 5 
and 10 pc from the star cluster centre, respectively, both at 
tsN = 5 x 10 5 yr. Figure [t] shows the soft X-ray luminosities 
for these two models compared to a model with a SN placed 
at the cluster centre (sr0tsn5zv). 

We can see that the X-ray luminosities increase when 
the SN explodes closer to the edge of the bubble. The max¬ 
imum soft X-ray luminosity increases by a factor of ~ 1.25 
between the model with supernova explosion at f?sN = 0 
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0 2x10 5 4X1 Cr 6x1 o 5 8x1.0 s 
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Figure 5. Evolution of the total soft X-ray luminosities for 
models with the supernova explosion occuring at the centre of 
the star distribution. The red lines are the models with homo¬ 
geneous metallicity (for a supernova event at t=2, 3, 5, and 
7.5x10° yr dash-dotted, solid, dot-dasli-dotted and dashed lines 
respectively). The blue lines are the models with different metal- 
licities (Z IS m = 0.3 Zq, Z wind = 3 Z 0 and Z SN = 10 Zq for 
the interestellar medium, stellar wind and supernova explosion, 
respectively, for SN explosions at t=2, 3 and 5x10 s yr). 


Table 2. Maximum soft X-ray luminosity and time intervals 
in which the soft X-ray emission remains above 10 34 , 10 3s and 
10 36 erg s _1 . 


Model 

Lmax,soft 

At a , 34 

At a ,35 

Ais,36 


[erg s” 1 ] 

[10 4 yr] 

[10 4 yr] 

[10 4 yr] 

sr0tsn2z0.3 

2.98X10 35 

6.42 

3.55 

— 

sr0tsn3z0.3 

2.83x10 s5 

5.60 

3.05 

— 

sr0tsn5z0.3 

2.98x10 s5 

6.01 

2.84 

— 

sr0tsn2zv 

3.01x10 s5 

6.72 

3.63 

— 

sr0tsn3zv 

2.92x10 s5 

5.96 

3.13 

— 

sr0tsn5zv 

3.08x10 s5 

6.34 

2.92 

— 

sr0tsn7.5zv 

2.68x10 s5 

6.13 

2.90 

— 

sr0tsn2zvC 

3.65x10 s ® 

8.20 

4.30 

— 

sr0tsn5zvC 

4.03x10 s5 

9.00 

3.96 

— 

sr5tsn5zv 

4.16x10 s5 

7.61 

3.79 

— 

srl0tsn5zv 

1.23x10 s6 

13.81 

7.27 

3.02 


and Rsn = 5 pc and a factor of ~ 4 when the SN explodes 
near the cluster radius (Rc = 10 pc), reaching a maximum 
luminosity of L max = 1.23 x 10 s6 erg s _1 . 

For model srl0tsn5zv, the only one that reached 
10 36 erg s“\ the time interval spent above 10 36 erg s _1 was 
30 kyr. In addition, this last model predicts a time spent 
above 10 35 erg s _1 of 72 kyr, and one above 10 34 erg s -2 
of 14 kyr. These numbers are ~ 3 times larger than those 
of the model with the supernova explosion occurring at the 
centre of the star cluster. 

From these results we can see that, on the one hand, 
for the case of the SN at the cluster centre the soft X-ray 
luminosity increase is not very sensitive to the time at which 
it is detonated. The luminosity increase is only slightly larger 


R sn =0 [pc] (w/cond) 



Figure 6. Evolution of the total soft X-rays luminosities for mod¬ 
els with supernova explosion in the centre of the star distribution 
with inhomogeneous metallicity. The magenta line are the models 
with thermal conduction and the blue lines are the models with¬ 
out thermal conduction process (for supernova event at t=2, and 
5 X 10 s yr dashed and solid lines, respectively). 


t s «=5x 1 0 5 [yr] 



Figure 7. Evolution of the total soft X-rays luminosities for mod¬ 
els with supernova explosion at t=5xl0 5 yr, and with inhomoge¬ 
neous metallicity. The blue,red and olive lines are the model with 
supernova event in R=0, 5 and 10 pc, respectively. 

for SN explosions that occur later in the evolution of the 
superbubble. On the other hand, when the SN is off-centre 
the luminosity increase depends considerably on the distance 
to the centre of the cluster. 

The color maps of Figure[8]show the soft X-ray emissiv- 
ity in the same layout as in Figure[3] Note that the blowout 
region (located at the upper right corner of the simulation 
box) provides the largest contribution to the luminosity in¬ 
crease seen 30 kyr after the explosion (middle panels). By 
200 kyr after the explosion (bottom panels), the emissiv- 
ity has almost returned to values comparable to its pre-SN 
level; however, due to the now larger X-ray emitting volume, 
the luminosity remains slightly above the original level (see 
Figure [7]). 

An interesting exercise is to compare the predicted lu¬ 
minosity statistics of our simulations with those of observed 
superbubbles. For this purpose we have taken a sample of 
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e)sr0tsn5zvC 
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f)sr10tsn5zv 
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Figure 8. Comparison of the X-ray emissivity in the soft (color maps, with the shown logarithmic scale given in erg s —1 cm -2 sterad" 1 ) 
and hard (contours, logarithmically-spaced levels from 10’ 11 to 10“ 8 erg s -1 cm -2 sterad -1 ) bands at the three different stages shown 
in Figure [3] 


26 bubbles with luminosities greater than 10 34 erg s 1 from 


the literature (OeyJ 1996 

Jaskot et al. 2011 Reyes-Iturbide 

et al. 2014 Dunne et al. 

• TZ TTZZ _ ' 

2001). Out of these, 18 have lu- 

T=T „ i__ o _ ^1_ j 


with L so ft > 10 3<! erg s _1 . We can use our numerical re¬ 
sults to predict how many bubbles out of these 26 should 
have luminosities above these two levels. In order to do this, 
we computed, from Table 2, the ratios of the times spent 
above these levels to the overall time where luminosity is 
above 10 34 erg s _1 for model srl0tsn5zv (the only one that 
reaches 10 36 erg s _1 ). We find that A t s ,36 / At s ,34 ~ 52.6 % 
and At s , 35 /Af s ,34 ~ 21.9 %. Thus, assuming that all bubbles 


in this observed sample reach a luminosity of 10 36 erg s ~ 4 
at some point in their evolution, this model predicts that 
about 14 should have a luminosity above 10 35 erg s -1 while 
about 6 should be observed above 10 36 erg s _1 . Though the 
values do not coincide exactly, they reasonably match the 
luminosity ratios of L 36 /L 34 ~ 12 % and L 36 /L 35 ~ 70 % 
observed in the sample. Here, L n is a soft X-ray luminosity 
that is greater than or equal to 10 " erg s _1 . 

There could be several explanations for this difference. 
For one, it is hard to judge whether this small sample of su¬ 
perbubbles is representative of the general population, and 
thus some variability can be expected in the statistics. At the 
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Table 3. Maximum hard X-ray luminosity and time inter¬ 
vals in which the hard X-ray emission remains above 10 33 and 
10 34 erg s _1 . 


Model 

Lmax,hard 

[erg s -1 ] 

Ath,33 

[10 4 yr] 

Ath,34 

[10 4 yr] 

sr0tsn2z0.3 

2 . 88 X 10 34 

2.00 

1.35 

sr0tsn3z0.3 

2 . 88 X 10 34 

1.99 

1.35 

sr0tsn5z0.3 

2 . 88 X 10 34 

1.96 

1.33 

sr0tsn2zv 

1.02X10 35 

2.03 

1.84 

sr0tsn3zv 

1.02X10 35 

2.02 

1.85 

sr0tsn5zv 

1.02X10 35 

2.00 

1.83 

sr0tsn7.5zv 

2.81X10 32 

0.00 

0.00 

sr0tsn2zvC 

9.90X10 34 

2.03 

1.84 

sr0tsn5zvC 

9.93X10 34 

2.00 

1.82 

sr5tsn5zv 

9.07xl0 34 

2.00 

1.81 

srl0tsn5zv 

9.84X10 34 

3.57 

1.99 


same time, our numerical results suggest that the position 
of a supernova explosion occurring inside the bubbles deter¬ 
mines whether a luminosity of 10 36 erg s _1 is reached at all 
during their lifetimes. Thus, if not all of the observed bub¬ 
bles have had off-centre SN explosion, it is to be expected 
that fewer of them would be observed above 10 36 erg s _1 
than what our models predict. 

3.2 Hard X-ray emission 

Hard X-ray emission is produced in the hottest regions in¬ 
side the bubble where individual winds interact, when the 
gas flow is faster than a ~ 1000 km s , as well as during 
the early stages of the SN remnant evolution and where the 
cluster wind and/or the SN remnant are heated by the re¬ 
verse shock. Thus one should expect that metallicity should 
have a significant effect on the hard X-ray emission. 

In Table [3] we show the maximum thermal hard X-ray 
luminosity, and the time intervals for which the luminos¬ 
ity remains above 10 33 erg s^ 1 (Ath, 33 ) and 10 34 erg s _1 
(Ath,33)- 

Figure [9] shows the evolution of the thermal hard X-rays 
luminosity for all models with the SN placed at the centre. 
In the models with different metallicities all three compo¬ 
nents exhibit maximum luminosities (~ 3 x 10 34 erg s _1 ) 
that are ~ 3 times larger than those of the models with 
homogeneous metallicity (~ 1.0 x 10 34 erg s _1 ). We have 
to mention that the thermal hard X-ray emission produced 
inside the star cluster (regions of wind collisions) is under¬ 
estimated by the models with homogeneous metallicity due 
to the rather low metallicity of the wind sources (0.3 Zq). 
In the models with the variable metallicity the wind is in¬ 
jected with a more appropriate metal content (3 Zq), thus 
the hard X-ray luminosity in these models should be closer 
to reality. 

We can see from Table [3] that the maximum hard X-ray 
luminosity is significantly larger in the models with variable 
metallicity, typically an increase of ~ 3.5 times. The time 
interval that the emission remains above 10 33 erg s _1 is simi¬ 
lar, although slightly larger in the models with more realistic 



Figure 9. Same as Figure [ 5 ] for hard X-rays luminosities. 


R sn =0 [pc] (w/cond) 



Figure 10. Same as Figure[6]for hard X-rays luminosities. 


metallicity. In contrast, the time that the emission remains 
above 10 34 erg s -1 is much larger (a factor of ~ 4) than that 
obtained in the models with homogeneous metallicity. 

An important fact to notice is that the ratio of the max¬ 
imum luminosities in soft-Xrays to hard X-rays is on the or¬ 
der of 10. Velazquez et al. (2013) explored the ratios between 
soft and hard X-ray emission in the young star cluster M 17. 
This particular cluster is partially immersed in the cluster 
parental cloud, and their models resulted in a ratio of soft 
to hard X-rays of two orders of magnitude. From our models 
we notice that the time intervals of high luminosity for soft 
X-rays (Af S[ 34 and Af a! 35 ) are larger than those obtain for 
hard X-rays (Ath ,33 and Ath, 34 ) 

This shows that very young star clusters with a SN 
event can produce hard X-ray luminosities that are only an 
order of magnitude fainter than the soft X-rays. However, 
this happens only for a short interval at the earlier stages of 
the SN remnant. After that the hard X-ray emission drops 
abruptly. 

Figure [lO] shows the evolution of the hard X-ray lumi¬ 
nosity for the two models with thermal conduction and their 
counterpart without thermal conduction. The maximum lu¬ 
minosities and the time intervals of high hard X-ray emission 
are remarkably similar in spite of the thermal conduction. 
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t slJ =5x 1 0 5 [yr] 



Figure 11. Same as Figure [7] for hard X-rays luminosities. 


Although small differences can be seen, the position of the 
supernova explosion in the cluster does not have a significant 
influence in the overall and maximum thermal hard X-ray 
luminosity (see Figure 111. 

It is also intresting to note that for an SN exploding at 
the centre of the cluster the highest luminosity achieved is 
much lower if the explosion occurs at later times (c.f. the 
long dashed curve in Figure [9] to the others). 

The distribution of the hard X-ray emission before the 
SN explosion, at peak luminosity, and after the luminosity 
returns to its previous level can be seen as the (logarithmi¬ 
cally spaced) contour levels in Figure[8] The emission during 
peak luminosity (middle row) is slightly more extended in 
the case where the SN is detonated off-centre (panel f), but 
returns to being concentrated inside the star cluster after 
the effect of the explosion has had time to decay (bottom 
row). The impact of thermal conduction on the hard X-ray 
emission is also evident. As can be seen, before (top row) or 
some time after (bottom row) the SN explosion the hard- 
band emission is at a much higher level in the cases with¬ 
out thermal conduction (left and right columns). However, 
during the luminosity peak (middle row) all models display 
important hard X-ray emission regardless of the inclusion of 
thermal conduction. The difference lies mainly in that in the 
case including thermal conduction (panel e) the emission is 
slightly more centralised than in the cases without. 

We must note that the resolution used in the models is 
not enough to capture all the details of the flow. The use of 
higher resolution allows larger compression factors as well 
as more small scale structure inside the bubbles. To esti¬ 
mate the uncertainty in the X-ray luminosities due to poor 
resolution we have taken a test case (model sr0tsn5z0), and 
reproduce the setup in the Walicxe-3d code ((Toledo-Roy et| 
ah] 2014). This code has adaptive mesh refinement (AMR), 


which allows to increase the resolution at a lower computa¬ 
tional cost, but the thermal conduction is not fully imple¬ 
mented. We ran the test case at an equivalent resolution of 
512 3 and 1024 3 cells, and while the details of the flow are dif¬ 
ferent, the integrated X-ray luminosities seem to reach con¬ 
vergence. The peak luminosity in the higher resolution runs 
is a factor of ~ 2 larger than in the 256 3 model sr0tsn5z0. 
And the times above 10 34 , 10 35 , and 10 36 erg s _1 are larger 
by a factor of ~ 1.5. All the results presented above have 


an uncertainty of this order of magnitude due to the limited 
resolution. 


4 COMPARISONS WITH OBSERVATIONS 


It is useful to compare our numerical models to the observa¬ 
tions of particular superbubbles. We have thus turned our 
attention to four superbubbles located in the Large Magel¬ 
lanic Cloud (LMC): N70, N185, DEM L50 and DEM L152. 
These superbubbles have some particular features (as we 
will show) that make them compatible with our results. 

N70 is a superbubble with a radius of approximately 
53 pc. According to the observations there is a SN located 
closer to the center than the edge of the cluster. The X- 
ray luminosity reported by Reyes-Iturbide et al. (20141 is 
about 2.4(±0.4) x 10 35 erg s _1 . We can see that our numer¬ 
ical results are in good agreement with these observations, 
in particular the models with the SN explosion in the center 
of the cluster. 

The case of N185 is quite similar to that of N70. N185 
has a spherical shape with an approximate radius of 43 
pc (Oey 1996). The X-ray luminosity obtained by Reyes- 

C 1 _ _ . 71 . . oe i ^-*- 


Iturbide et al. (2014) is 2.1 (±0.7) x 10 35 erg s . Following 


Rosado et al. (19821, a possibility to explain the high veloc¬ 


ity of this superbubble is that a SN explosion ocurred. From 
its spherical shape, we conclude that the SN explosion must 
be located near the centre of the cluster. As in the previous 
case, our numerical models are in good agreement with the 
observed X-ray luminosities. 

Two other interesting cases are DEM L50 and DEM 
L152. These are two superbubbles with very intense X-ray 
emission. According to observations DM L152 has a radius 
of approximately 50 pc ( Jaskot et al.j|2011 l, and DEM L50 
has roughly the same radius (Oey 1996). Jaskot et al. (20111 
reported an X-ray luminosity in the range 2.0-4.0xl0 JO erg 
s _1 for DEM L50 and an emission of 5.4-5.7xl0 35 erg s _1 . 
These two superbubbles contrast with N70 and N185 in that 
they contain an off-centre SN explosion, which is clearly 
distinguishable in the observations. Only numerical model 
snl0tsn5zv predicts a luminosity comparable to the observed 
values. The luminosities predicted by other models are not 
high enough to match with these observed values. 

In order to better compare with observations, we have 
calculated the hardness ratio for three of our models (see 
Figure 121. These models are those discussed in Section 3: 
two with a centred SN explosion that differ in the presence 
of thermal conduction, and the third with the SN explosion 
10 pc off-centre. Following Jaskot et al. (2011), the hardness 
ratio is defined as: 

H-S 


Q = 


H + S 


( 12 ) 


where H is the flux energy in the 2-10 keV energy band 
(corresponding to hard X-rays) and S is the flux energy in 
the 0.1-2 keV energy band (corresponding to soft X-rays). 
After computing the fluxes and obtaining Q for each cell 
in the simulation, we integrate along the z axis in order 
to project the result into a 2D map, assuming that the X- 
ray absorption due to the material inside the bubble can be 
neglected. 

In Figure [l2| we observe that, before the SN explosion 
occurs, in models without thermal conduction (left and right 
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Figure 12. Hardness maps for the same cases shown in Figure [3] 


columns) the hardness ratio peaks at ~ —0.55 at the cen¬ 
tre of the stellar distribution, and decreases as we approach 
the edge of the cluster (the shell of swept-up ISM material 
emits mainly soft X-rays). In the case of the model that 
includes thermal conduction (middle column), we observe 
that Q ~ —1. in most of the bubble, indicating that hard 
X-ray emission is largely negligible. As a result, the mate¬ 
rial ejected by the stars is cooled quickly from hard X-ray 
emitting temperatures (10 8 K) to temperatures in the range 
10 5 - 10° K where soft X-ray emission is favored. 

The SN explosion drastically alters the hardness maps 
for the models without thermal conduction. The SN shock- 
wave sweeps the cluster volume, devoiding the center of the 
bubble of hard X-ray emitting gas and forming regions with 


Q ~ —0.6 closer to the edge of the bubble. In the model 
with thermal conduction, the hard X-ray emission is small 
to begin with, and the effect of the explosion on the hardness 
ratio is not as noticeable. 


The predicted Q values obtained in our numerical mod¬ 
els are similar to those obtained by Jaskot et al. (20111 for 
DEM L50 and DEM L152. Nevertheless, the specific details 
and assumptions of our simulations make it hard to establish 
direct comparisons to specific observed bubbles. In order to 
use the hardness ratio to predict some of the physical pro¬ 
cesses that occur or have occurred in the super-bubbles we 
would need to separately simulate the specific physical de¬ 
tails, such as the position and mass and energy injection 
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rates of each star, the ISM density, of each particular bub¬ 
ble, which is out of the scope of this work. 


5 CONCLUSIONS 


In this paper we present 3D hydrodynamical models of 
the evolution of the soft and hard thermal X-ray luminosi¬ 
ties produced inside superbubbles driven by massive stellar 
winds including the effect of a supernova explosion and ther¬ 
mal conduction. 

In all models we include the injection of mass and en¬ 
ergy by a cluster of wind sources and a single SN event. We 
have varied the position of the SN with respect to the centre 
of the cluster as well as the detonation time. Also, we have 
worked out models with a uniform metallicity and models 
in which the environment, the winds and the SN have dif¬ 
ferent metallicities. The metallicities are used to calculate 
the radiative cooling rate, and has an effect on the emissiv- 
ity in X-rays. We have also taken into account the effects of 
thermal conduction in two of the models. 

In the models with different metallicities we used Z = 
0.3 Zq for the environment, 3 Z@ for the winds and 10 Zg 
for the SN ejecta. Our models show that the contribution of 
the metallicity of the winds and the supernova remnant is 
negligible for the soft X-ray emission of superbubbles, but 
becomes important for the hard X-ray component. In these 
models the ratio of soft to hard maximum luminosity can 
be as extreme as 10 (i.e. the hard X-ray luminosity reaches 
10% of the one for soft X-rays). 

The models with thermal conduction result in a notice¬ 
able increase in the total luminosity of soft X-rays, by a 
factor of ~ 1.25. However, this factor is smaller than the 
two orders of magnitude difference predicted in the stan¬ 


dard model of Weaver et al. (1977 I and Chu & Mac Low 


(1990). The differences are likely to come from the fact that 


the standard model of Weaver et al. (1977) considers just 


a single star, and the extension to a star cluster in |Chu fe| 
Mac Low (19901 and Chu et al. (19951 does not account for 


the cooling of the gas due to the interaction of the stellar 
winds. Thermal conduction has a slightly larger effect on 
the total integrated emission of hard X-rays, increasing the 
luminosity by a factor of ~ 2.6. 

The most important contribution to the emission of soft 
and hard X-rays is produced by the injection of mass and 
energy by supernova explosions. In soft X-rays the luminos¬ 
ity increases by up to two orders of magnitude when we 
consider a supernova explosion placed at the cluster centre, 
and up to three orders when it explodes at the edge the star 
cluster. 

Another important factor to consider is the time dur¬ 
ing which the luminosity remains high (i.e. observable). We 
show that when off-centre supernova events occur (close to 
the shell) the luminosity can increase by one or two orders 
of magnitude above that predicted by the standard model 
without SN, and that it can be maintained by a few tens 
of thousands of years. Indeed, as the supernova explosion 
occurs closer to the shell of swept up ISM, the maximum 
luminosity of soft X-rays as well as the time interval during 
which luminosity is enhanced increase. 

An important increase in the maximum soft X-ray lu¬ 
minosity is produced when the SN ejecta collide with the 


dense shell of swept up ISM gas left behind by its interac¬ 
tion with the cluster wind. In these cases X-ray luminosities 
of 10 46 erg s -1 can be achieved. On the other hand, su¬ 
perbubbles where the SN explosions have not taken place 


near the shell, such as N 70 and N 185 (Jansen et al. 2011 
and Reyes-Iturbide et al.||2014 ), have lower X-ray luminos¬ 
ity, and can be explained using our models with a slightly 
off-centre SN. 

In clusters without SN events, or with a SN placed at 
the centre of the cluster, the contribution to the luminosity 
made by the SN is hard to observe, in particular because 
the observable flux increase in the soft X-ray emission lasts 
for a short time. This could be happening in massive stel¬ 
lar clusters in the Galaxy, such as Arches, Quintuplet and 
NGC3603, that have a hundred massive stars with a total 
observed X-ray emission of ~ 10 44 erg s _1 . 
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